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Abstract 



The ground state rotational spectrum of PH3 has been reanalyzed taking into account recently published very accurate data from 
C*~) sub-Doppler and conventional absorption spectroscopy measurements as well as previous data from the radio-frequency to the far- 
infrared regions. These data include A J = AK = transitions between A\ and A 2 levels, A J = 0, AK = 3 transitions as well as 
regular A J — 1, AK — rotational transitions. Hyperfine splitting caused by the 31 P and 'H nuclei has been considered, and the 
treatment of the A1/A2 splitting has been discussed briefly. Improved spectroscopic parameters have been obtained. Interestingly, 



the most pronounced effects occured for the hyperfine parameters. 
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Phosphine, PH3, was detected as a trace constituent of the 
atmospheres of the giant planets Jupiter [1] and Saturn fl, re- 
spectively by infrared spectroscopy. Observations of the lowest 
7=1-0 rotational transition near 267 GHz toward Jupiter 
and Saturn [4] were also reported. Using a far- and mid-infrared 
(FIR and MIR) spectrometer on board of the Cassini satel- 
lite, it was even possible to study the meridional PH3 distri- 
bution in Saturn's upper troposphere fl. Higher resolution FIR 
observations of Saturn's atmosphere were carried out very re- 
cently with the Herschel satellite, and the vertical distribution 
of PH3 was derived |@]. Phosphine is also a likely species 
to occur in the atmospheres of giant extrasolar planets or of 
brown dwarfs [7]. Moreover, it has been tentatively detected in 
the circumstellar shells of carbon-rich late-type stars CW Leo- 
nis lHH, also known as IRC +10216, and V1610 Cygni, also 
known as CRL 2688, by radio-astronomical observations of the 
7=1-0 transition. 

Small molecules, and even more so light hydride species, 
such as PH3, are also of great interest for basic science. In the 
early years, an important questions was if PH3 would display 
measureable inversion splitting as its lighter homolog ammo- 
nia, NH 3 , does E[H. 

Early reports on the ground state rotational spectrum of PH3 
were published by the Gordy group 11121 1 1 311 . Burrus fl 1 411 de- 
termined the dipole moment from Stark measurements. Davies 
et al. Il5ll studied transitions between A\ and A2 levels with 
AJ = AK — at radio-frequencies (RF). They were unable to 
resolve any tunneling splitting, but obtained hyperfine structure 
parameters caused by the 31 P and 'H nuclei and an improved 
value for the dipole moment. Chu and Oka recorded AK = 3 



transitions with AJ = and K = +1 - +2 and estimated an 
effective value for the centrifugal distortion induced dipole mo- 
ment from intensity measurements Ill6fl. Later, corresponding 
transitions were reported for K = 0-3 [17], K = 2-5 [18], and 
K — 3 - 6 [ 19]. Regular rotational transitions with AJ — 1 and 
AK — in the millimeter and sub millimeter regions M 1 8L Il9 1 



'Corresponding author. 

Email address: hspmOphl . uni-koeln . de (Holger S.P. Miiller) 



and in the far-infrared (FIR) region [20] were also studied. 

Recently, Cazzoli and Puzzarini 112 111 obtained very accurate 
transition frequencies for the 7=1-0 and 2-1 transitions em- 
ploying sub-Doppler (Lamb-dip) spectroscopy, which resolved 
31 P and 'H hyperfine structure, and for the 7 = 3-2 transi- 
tions using conventional absorption spectroscopy. As this study 
employed in their fit very limited data beyond their own, a 
combined fit of these data together with previous results has 
been carried out here to obtain an improved set of spectro- 
scopic parameters which should be useful not only for spec- 
troscopic observations of phosphine in planetary atmospheres, 
brown dwarfs, or other regions in space, but also for future anal- 
yses of rovibrational spectra or as benchmarks for quantum- 
chemical calculations. 

2. Spectral analysis and discussion 

Phosphine is an oblate rotor with a permanent dipole moment 
of 0.57395 (3) D lfl5ll . It is moderately close to the spherical 
rotor limit with y = (Co - Bq)IAq = -0.1361 being small in 
magnitude, but not particularly so. 

The rotational energies and transitions of a strongly oblate 
rotor in its ground vibrational state and without Ai/Aj splitting 
require for their modeling only C - B and B along with distor- 
tion parameters, Dj, D jk, and Dk to lowest order. In a strongly 
prolate rotor, C is replaced with A. The purely axial parameters 
C - B, Dk, etc. can usually not be determined by rotational 
spectroscopy. Their determination requires, e.g., the observa- 
tion of centrifugal distortion induced AK = 3 transitions. The 
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Table 1: Spectroscopic parameters" (MHz) of phosphine, PH3, from the present fits in comparison to selected previous fits. 





Present, preferred 6 


Present, alternative^ 


Fusina and Carlotti'' 


Cazzoli and Puzzarini e 


Davies et a\J 


B 


133480.12821 (15) 


133480.127 94 (15) 


133480.122 80 (177) 


133 480.128 989 (95) 


— 


C 


117 489.4304 (73) s 


117489.427 (73) s 


117489.4256 (73) 


117 489.436 s -' 1 


— 


Dj 


3.936391 (22) 


3.936210 (22) 


3.936059 (105) 


3.936901 (36) 




Dnr 


-5.168 535 (78) 


-5.167 110 (78) 


-5.167 280 (86) 


-5.171 02 (44) 




D K 


4.23445 (54) 


4.232 60 (54) 


4.233 82 (51) 






Hj x 10 s 


348.906(142) 


326.315 (142) 


325.74 (39) 


415.7 (39) 




H JK x 10 6 


-868.02 (50) 


-664.97 (50) 


-665.33 (73) 


-1 237. (85) 




H KJ x 10 6 


801.49(148) 


463.43 (149) 


464.10(276) 


1 340. (270) 




H K x 10 6 


-31.4(111) 


114.8 (111) 


152.7 (103) 






hi x 10 6 


11.3* 


_C 








Lj x 10 9 


-37.667 (218) 


-37.672 (218) 


-36.96 (43) 






Ljjk x 10" 


113.71 (105) 


149.91 (106) 








L JK x 10 9 


-222.0 (33) 


-408.7 (33) 


349.0 (39) 






Lkkj x 10' 


247.9 (50) 


519.5 (50) 


-547.7 (93) 






ex 10 3 


714.508 2(42)' 


831.6522(36)' 


829.254 (56) 






£j x 10 6 


-168.896(72)' 


-180.210(62)' 


-179.611 (204) 






e K X 10 s 






265.9 (61) 






Cjj x 10' 


37.25 (29)'' 


34.60 (24)'' 


34.97 (199) 






e JK x 10' 






-71.5(102) 






€ K K X 10' 






87.2 (82) 






C N (P) x 10 3 


114.897 (15) 


114.897 (15) 




115.35 (12) 


114.90(13) 


C K (P) x 10 3 


116.396 (35) g 


116.396 (35) g 




115.0 (14) s 


116.38 (32) 


C N (H) x 10 3 


-7.998 (9) 


-7.998 (9) 




-7.57 (13) 


-8.01 (8) 


C K (H) x 10 3 


-7.704 (22)f 


-7.704 (22)f 




-7.69 s '' 


-7.69 (19) 


-3Di(P-H)x 10 3 


3.01 (35) 


3.01 (35) 




3.05 (62) 


k 


-0.5D 2 (P-H) x 10 3 


-9.15' 


-9.15' 




-9.15' 


k 


1.5£> 3 (H-H)x 10 3 


20.62(13) 


20.62(13) 




25.3 (19) 


k 



°Numbers in parentheses are one standard deviation in units of the least significant figures. 

*The parameter /13 was kept fixed to the value calculated in Ref. [ 231 . 

c The parameter A3 was constrained to zero. 

rf Ref. [ 20], preferred fit with /13 etc. kept fixed to zero. 

e Ref. H. 

^Ref. Q. 

s The parameters C was determined from C - B and B, taking correlation into account. Analogously, Ck was calculated from Ck - Cn and Cn . 

h C was kept fixed to the value in Ref. [22]. 

'The parameters determined in the fit are V2e etc. 

^Ck(H) was kept fixed to the value in Ref. H3l . 

*The tensorial spin-spin coupling parameters were kept fixed to values calculated from structural parameters, but were not given specifically. 
'£>2 was kept fixed to the value calculated in Ref. |21]. 
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A1IA2 splitting is modeled to lowest order by h^J^ + J 6 } with 
J ± = J x + iJ y . J and K distortion corrections and possibly even 
higher order terms may be required in the fit. Early treatments 
of the rotational spectrum of PH3 used parameters as described 
above to fit the rotational transitions. However, it was recog- 
nized early on that the parameter derived for I13 was too large 
by about a factor of 4 II 1 511 . Moreover, it was found later, that 



the Hamiltonian converged rather slowly [19]. 

j — 1 

Tarrago and Dang Nhu [22] discussed the formulation of 
the Hamiltonian for symmetric top molecules approaching the 
spheric rotor limit employing the examples OPF3 and PH3. 
They pointed out that for such molecules not only the AK = 6 
terms derived from h-}, also called in their work, affect the 
A1/A2 splitting, but also terms with AK = 3, whose lowest or- 
der term is e, which is also called ^3. The parameter e is the 
coefficient of [/+ + Jl , J z ]+, were [A, B]+ = AB + BA is the anti- 
commutator of A and B. They presented one fit with parameters 
derived from /z3, and one fit with parameters derived from e and 
h?,. The latter fit was slightly better, but yielded a value for /13 
only slightly smaller than in the fit without e, meaning that then- 
value was still too lage by almost a factor of 4. Moreover, then- 
value for e was too small by about a factor of 3. Interestingly, 
they did not attempt any fits without parameters derived from 
h 3 . 

Fusina and Carlotti recorded the FIR spectrum of PH3 and 
made assignments for the ground vibrational state between 44 
and 202 cirT 1 with J" between 4 and 22 and K between 19 and 
ll20ll . They used in their fits their o wn data together with previ- 
ous data from Refs. |fl3j [Hj fl6[ |l7l 
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1911 . Three different fits 
were presented which included parameters derived from both e 
and /13 as well as parameters derived either from e or /13. The 
parameters e and I13 were completely correlated if both were 
used. The parameters /13 and h^K were completely correlated 
in the fit without parameters derived from e. In contrast, the 
fit with only parameters derived from e did not produce severe 
correlations, yielded the best standard deviation, but employed 
two more parameters than the other two fits. This fit was their 
preferred one. 

The fits carried out in the course of the present investiga- 
tion employed Herb Pickett's spfit and spcat programs [23]. 
They started from the preferred parameters from Fusina and 
Carlotti 1I20I1 and the hyperfine parameters employed by Cazzoli 
and Puzzarini lpH] using transition frequencies from Refs. lfl5l 
[l ^ fl ^ fl^lT ^l iallll] with uncertainties as reported. The RF 
transitions 01511 were reportedly accurate to between 0.1 and 
0.2 kHz. Initially, 0.2 kHz were used in the fits because the 
published deviations from the measured frequencies exceeded 
0.6 kHz in two cases. However, the transitions were reproduced 
in the present fits to well within 0. 1 kHz on average with no se- 
vere outliers. Therefore, 0.1 kHz uncertainties were employed 
for these data in the final fits. The uncertainties of Ref. [18] 
were used as reported in Ref. IU9II . Including Ljjk in the fit, 
which was contrained to zero in Ref. lEoll . yielded a reasonable 
value, and suggested the parameters ex, £jk, and €kk to be in- 
determinate. The latter three parameters could be omitted from 
the fit without significant deterioration of the fit. The result- 
ing parameters, including those for the hyperfine structure, are 



given in Table Q]under heading "Present, alternative". 

Hyperfine structure caused by the 31 P and 'H nuclei was re- 
solved in the RF spectra lfl5ll and in the Lamb-dip measure- 
ments from Ref. 112 ill . The hyperfine Hamiltonian was de- 
scribed in sufficient detail e.g. in the latter work. As mentioned 
in that work, the Lamb-dip measurements were rather sensitive 
to the perpendicular components Cm of the P and H nuclear 
spin-rotation tensors, which are also known as Cj_, and to the 
axial components Dj(P-H) and £>3(H-H) of the spin-spin cou- 
pling tensors, which are also known as ii||(P-H) and (H— H), 
respectively. They were not as sensitive to the axial components 
Ck of the nuclear spin-rotation tensor, which are also known as 
C||, and to the perpendicular component £>2(P-H) of the P-H 
spin-spin coupling tensor, also known as d ± (P-H). In contrast, 
the RF data were sensitive to Cn and Ck of both nuclei, but not 
as sensitive to the spin-spin coupling tensors. In the present fit, 
D2 could still not be determined significantly and with confi- 
dence. However, its omission increased the standard deviation 
of the fit somewhat. Therefore, it was retained in the fit kept 
fixed to the calculated value from Ref. [21]. Other parameters, 
such as C_(H), the off-diagonal H spin-rotation parameter, the 
scalar spin-spin coupling parameters, or distortion corrections 
to the P or H nuclear spin-rotation parameters were not deter- 
mined with significance, and their inclusion had a negligible 
effect on the quality of the fit. Therefore, they were omitted 
from the final fits. 

As can be seen in Table Q] the present hyperfine parameters 
are in good to very good agreement with previously determined 
ones I T51 21 ], but have much smaller uncertainties in the present 
fits. The values of the spin-spin coupling parameters -3D; and 
1 .5D3 agree very well with the calculated ones of 3 . 1 8 kHz and 
20.7 kHz, respectively. In fact, keeping the spin-spin coupling 
parameters fixed to the calculated values, does not deteriorate 
the quality of the fit significantly. In addition, the values of the 
other parameters and, somewhat surprisingly, even their uncer- 
tainties are affected negligibly. 

The remaining spectroscopic parameters agree for the most 
part well with the preferred ones from Fusina and Carlotti lEoll 
as far as they were determined in both fits. Exceptions are the 
two octic distortion parameters Ljk and Lkkj, which is a re- 
sult of the omission of Ljjk by Fusina and Carlotti and of the 
omission of some higher order distortion corrections to e in the 
present fits. Notably, the values of e agree quite well, about 
832 kHz versus 829 kHz. Both values, however, are slightly 
larger than 723 kHz calculated by Tarrago and D ang Nhu ||22 1 
and 773 kHz calculated by Fusina and Carlotti [20]. The pa- 
rameter e and its distortion corrections cause AK = 3 transi- 
tions, among others, to have non-zero intensities very similar 
to those calculated without these parameters, but with the ef- 
fective centrifugal distortion correction to the permanent dipole 
moment mentioned in Ref. fltSll . Therefore, another fit has been 
carried out with kept fixed at the value of 11.3 Hz calcu- 
lated in Ref. I22I1 . The e value of 714.5 kHz resulting from 
this fit is very close to the one calculated by Tarrago and Dang 
Nhu 1 22], but somewhat smaller than the one from Fusina and 
Carlotti ll20ll . All other parameters are also affected in their 
magnitudes. This fit is considered to be the preferred one. The 
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effects on the intensities of the AK = 3 transitions correspond 
closely to the square of the ratios of the e values, the intensities 
of all other transitions change only marginally between the two 
data sets. 

The quality of the preferred fit and of the alternative fit are 
identical. The experimental transition frequencies have been 
reproduced within the reported uncertainties on average, the 
rms error is 0.93, slightly better than the ideal 1.0. There 
is some scatter among the rms errors of individual data sets: 
it is 0.77 for the RF data |B and 0.64 (11, 0.87 and 



0.94 1181 11911 for the older millimeter and submillimeter data. 
The FIR data from Fusina and Carlotti [20] were given uncer- 
tainties of 0.0001 cirT 1 throughout, as suggested by the authors, 
even though the residuals increased somewhat for the transi- 
tions with the highest two J values. Two transitions omitted 
already by Fusina and Carlotti were also omitted here. In addi- 
tion, the frequency of the transition with J" = 16 and K = 9 
was omitted as its residual exceeded Act. The rms error of the 
FIR transitions was 0.94; the frequencies of 4 transitions, all 
at higher / values, had residuals larger than 3cr, which is rea- 
sonably close to the statistical expectations in the case of 230 
different transition frequencies. An rms error of 1.37 was ob- 
tained for the data from Ref . 112 111 . Even though this suggests 
slightly optimistic uncertainty estimates, no adjustments were 
made, as one frequency, 533815.1259 (10) MHz, contributed 
particularly to the rms error. Omission of this line reduced the 
rms error of this data set to slightly below 1.1. In addition, the 
omission of this line from the overall line list affected all pa- 
rameters within experimental uncertainties. Therefore, this line 
was retained in the final fit. 

The /-dependent parameters, in particular the ones depen- 
dent purely on J have much smaller uncertainties in the present 
fits compared to the preferred fit of Fusina and Carlotti 1I20I . 
This is hardly surprising given the nature of the new data from 
Ref. l2in . The uncertainties of Dk and Hk here are actually 
slightly larger than from the FIR fit. This may be caused by dif- 
ferences in the correlation of the parameters or by differences 
in the weighting scheme. 

Predictions of the rotational spectrum were made for the 
Cologne Database for Molecular Spectroscopy, CDMsQ [24, 
2511 . using the parameters of the preferred fit. The upper quan- 
tum numbers were J — 40 and K = 38, ensuring convergence of 
the rotational part of the partition function to 10 4 up to 500 K. 
The intensity cut-offs were chosen to include all experimental 
AJ = AK = transitions between A 1 and A2 levels on one hand 
and to limit the number of transitions with uncertainties equal 
to or larger than 1 GHz, as spcat displays such uncertainties as 
999.9999 MHz or 0.0334 cm -1 . This leads to transitions with 
upper J and K of 34 and 3 1 , respectively. The K = transition 
with / = 8 - 7 is the most intense transition at 300 K. At 500 K 
and 1000 K, the corresponding transition with / = 11-10 
and / = 14 - 14, respectively, is the most intense one. With 
phosphine being only a minor constituent in the atmospheres 
of Jupiter and Saturn, the current predictions should therefore 



contain sufficient ground state rotational transitions not only at 
300 K, but also, after appropriate conversion of the intensities, 
for 500 K and likely for 1000 K. However, the neglect of ro- 
tational transitions in excited vibrational states may cause non- 
negligible errors already around 300 K, as can be seen below. 

Intensities of the transitions with AK = 3 should be viewed 
with some caution as indicated above. Predictions of transition 
frequencies with AK = should be reliable as long as the pre- 
dicted uncertainties do not exceed 5 MHz. Predictions of tran- 
sition frequencies with AK + should be viewed with caution 
if the predicted uncertainties are larger than 0.5 MHz. In these 
cases, the actual frequencies may differ from the predictions by 
more than five times the predicted uncertainties. 

The partition function values given in the documentation of 
the CDMS entrj0do not include vibrational contributions; they 
are about 2.4 % at 300 K and amount to about 18.5 % at 500 K. 
Classically estimated vibrational contributions are intended to 
be provided through a link in the documentation file. It should 
be mentioned that the HITRAN entry takes into account vibra- 
tional contributions to the partition function. Moreover, the in- 
tensities at 300 K have been converted to those at 296 K, which 
is the default temperature in HITRAN. 

Pickett's spfit and spcat programs 02311 use the total parity. 
This parity is therefore also used in the ascii table version of the 
PH 3 entry in the CDMS. It differs both from the Ai/A 2 parity as 
well as from the A + /A_ parity. A levels with - parity belong to 
A2 for even K and to A 1 for odd K. A levels with - parity belong 
to A_ for J + K even and to A+ for J + K odd. Additional confu- 
sion may occur because the + sign has been omitted, and tran- 
sitions with rounded A1/A2 splitting of less than 0.1 kHz have 
been merged. These transitions are thus between two A levels 
without parity designation. Moreover, since only 2 characters 
are available for each quantum number, K values > 10 with neg- 
ative parity are indicated by lower case characters; e.g. a2, a5, 
a8, and bl are used for -12, -15, -18, and -21, respectively. 
The entry in HITRAN provides their common parity labels. It 
should be also mentioned that a process is under way to create a 
database version of the CDMS within the framework of the Vir- 
tual Atomic and Molecular Data Centre (VAMDCf|[28]. The 
official release of a test version is still pending because very 
many entries need to be imported properly, and in many cases 
necessary auxiliary files need to be created. One of the goals of 
VAMDC is providing easily interpretable quantum numbers. 

Fairly recently, dipole moment and potential energy surfaces 
were reported, and selected ground state rotational energies 
have been provided l26ll . The deviations between their ener- 
gies and the present ones are fairly small. They drop at J — 10 
from 0.9164 cirT 1 for K = to 0.4720 cirT 1 for K = 9. The de- 
viations are larger at J = 20, where they drop from 3.43 12 cm -1 
for K = to 2.9680 cirT 1 for K - 9. Because of the extent of 
the experimental data, these deviations are probably predomi- 
nantly due to insufficiencies in the calculated potential energy 
surface. 



'http://www.astro.uni-koeln.de/cdms/ 



2 http://www.astro.uni-koeln.de/cgi-bin/cdmsinfo?file=e034501.cat 
3 http://www.vamdc.eu/ 
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More recently, a very extensive and very accurate line list 
has been created for PH3 12711 . The authors found a remarkably 
good agreement between their transition frequencies and the 
ones in the first CDMS entry which had been created in October 
2008 and which has been recalculated in February of 2013 by 
lowering the intensity cut-offs to provide weaker transition fre- 
quencies which in part access higher quantum numbers. These 
predictions correspond essentially to the ones derived from the 
alternative fit in the present work. The authors mentioned that 
the deviations are below 0.003 cirT 1 for the strongest 100 tran- 



sitions 02711 . As these transitions involve J values between 3 
and 14 and K values from to 9, which are all covered by 
experimental data, the deviations are probably again predom- 
inantly due to insufficiencies in the calculated potential energy 
surface. Moreover, differences in the predicted frequencies gen- 
erated from the preferred fot or the alternative fit are very s mall 
for these transitions. Even more remarkable are the deviations 
at low quantum numbers. The Jk = lo - Oo prediction from 
Ref. [27] is only 2.8 MHz lower than the experimental fre- 
quency as well as the frequency calculated from either present 
fit. The deviations decrease first with increasing J, until the 
frequency of the 80 - 7o transition is actually higher than the 
one predicted from the preferred fit. Similarly, the 2±\ - 2 T 2 
transition frequency from Ref. ll27ll is only 3.6 MHz lower than 
the one predicted from the preferred fit. Again, the deviations 
decrease first for increasing / until the 10±i — 10 T 2 transition 
frequency is higher in frequency. The agreement is still rather 
good if all transitions from the previous CDMS predictions are 
compared with those from Ref. 112711 ; the authors give an rms 
value of 0.076 cm -1 . Because of the limitations in the experi- 
mental data as well as in the Hamiltonian model, it may well be 
that the predicted transition frequencies from the present study 
are eventually less reliable that those from Ref. fl27ll . 

The intensities from Ref. [27] agree very well with either pre- 
diction from the present fits for AK = transitions. The value 
for jj 2 from the former work is 3.9 % larger than the experimen- 
tal value from Ref. 11511: this deviation is almost canceled by the 
neglect of the vibrational contributions to the partition function 
at 300 K in the present work. However, there are larger devia- 
tions between the calculated intensities of AK — 3 transitions; 
e.g., the intensity of the 2 ± \ - 2 T 2 transition from Ref. fl27ll is 
a factor of 1.85 weaker than the intensity from predictions gen- 
erated from the preferred fit; it is even weaker by a factor of 
2.5 if predictions generated from the alternative fit are used. 
Possible explanations for these deviations are insufficiencies in 
the calculations from Ref. 12711 . insufficiencies in the preferred 
Hamiltonian model, e.g. €k had to be constrained to zero, or the 
neglect of distortion corrections to the electric dipole moment. 
It may well be that each of these explanations is responsible in 
part. 

3. Conclusion 

Improved spectroscopic parameters of PH3 have been deter- 
mined from recent very accurate rotational transition frequen- 
cies combined with previous data. Predictions for the ground 
state rotational spectrum of PH3 generated from the preferred fit 



and employing the permanent dipole moment of 0.57395 D [ 15 ] 
will be provided in the catalog section of the CDMS0 [2H2I 
and in the new edition of HITRAN0. Transition frequencies 
with assignments, uncertainties, and residuals between the ob- 
served frequencies and those calculated from the preferred pa- 
rameter set as well as the predictions in the CDMS format will 
be provided as supplementary material. Line, parameter, and 
fit files along with various auxiliary files will be provided in 
the archive section of the CDMS catalog. The old predictions 
together with available auxiliary files will be available in the 
archive section of the CDMSn. 
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